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In this paper, we investigate, phonon effects on the optical properties of a spherical quantum 
dot. For this purpose, we consider the interaction of a spherical quantum dot with classical and 
quantum fields while the exciton of quantum dot interacts with a solid state reservoir. We show that 
phonons strongly affect the Rabi oscillations and optical coherence on first picoseconds of dynamics. 
We consider the quantum statistics of emitted photons by quantum dot and we show that these 
photons are anti-bunched and obey the sub-Poissonian statistics. In addition, we examine the effects 
of detuning and interaction of quantum dot with the cavity mode on optical coherence of energy 
levels. The effects of detuning and interaction of quantum dot with cavity mode on optical coherence 
of energy levels are compared to the effects of its interaction with classical pulse. 

PACS numbers: 42.50.Ct, 42.50.Ar, 73.21. La, 63.20.kd 



I. INTRODUCTION 

The fundamental system in cavity quantum electro- 
dynamics (cavity-QED) is a two level atom interacting 
with a single-cavity mode Q-H- Recent developments 
in semiconductor nano-technology have shown that ex- 
citons in quantum dots (QDs) constitute an alternative 
two-level system for cavity-QED application Q. There 
are many similarities between the excitons in QDs and 
atomic systems, such as the discrete level structures 
which is subsequent of three-dimensional confinement of 
electrons. On the other hand, there are also important 
differences, for example coupling to phonons, carrier- 
carrier interaction and surface fluctuation. Coupling of 
electrons to phonons plays a major role in QDs. The cou- 
pling of phonons to the QD provides a basic dephasing 
mechanism and thus marks a lower limit for the deco- 
herence In self-assembled QDs it is indeed the 

clastic phonon scattering (pure dephasing) which dom- 
inates the loss of coherence on a picosecond time scale 
at temperatures below 100K [7f. The effects of electron- 
phonon interactions on strong exciton-photon coupling in 
cavity-QED has been considered @. It has been shown 
@ that the phonon-induced damping of Rabi oscillations 
in a QD is a non-monotonic function of the laser-field in- 
tensity that is increasing at low fields and decreasing at 
high fields. 

QDs are also promising candidates for efficient, deter- 
ministic single photon sources [10|-[ll|. Then the QDs 
are important sources of non-classical light. For this kind 
of application an understanding of the coherence prop- 
erties of its optical transitions is of great importance. 
Therefore, there are two processes in optical manipula- 
tion with semiconductor QDs: coherent control of the QD 
exciton state (l2j and measurement of quantum statistics 
of emitted light with QD [l3[ . A theoretical investigation 
of exciton dynamics and the possibility of generation of 
non-classical light has been considered without taking 
into account the phonon effects (T3 |. 

In this paper, we investigate the effects of electron- 



phonon interactions on optical coherence and quantum 
statistics of light emitted by a pulse driven QD interact- 
ing with a cavity mode. The photon statistics from a 
driven QD under the influence of the phonon environ- 
ment has been considered recently [l5j . On the other 
hand, influence of phonons on incoherent photon emis- 
sion of a QD in the presence of pulse excitation had been 
considered [16(. We use the most widely studied model 
for phonon effects in QDs which accounts two electronic 
levels coupled to a laser pulse and to non-interacting 
phonons |17| . As mentioned, phonon interaction provides 
a dephasing mechanism for optically induced coherence 
on a time scale (a few picosecond) much shorter than 
for radiative interaction and recombination [3. Due 
to the different correlation time for a phonon reservoir 
(few picosecond) and for a radiative reservoir (several 
ten nanosecond) we restrict our attention to the time 
scales which dephasing effects due to the phonon system 
play an important role (with the radiative reservoir we 
mean a reservoir for photon system. The mentioned time 
scale relates to the decay time of cavity photons). Then 
we do not consider any damping effect on cavity mode 
and spontaneous emission. In our consideration the only 
damping effect is related to phonons. 

The paper is organized as follows: In section II we de- 
scribe the model Hamiltonian and master equation that 
allows to calculate the evolution of populations and co- 
herence of the energy levels. In section III we present the 
exciton dynamics and its coherence while driven with a 
laser pulse. The photon statistics and exciton dynamics 
of pulse driven QD interacting with a cavity mode is pre- 
sented and discussed in section IV. Section V is devoted 
to a summery and conclusion. 



II. THEORETICAL MODEL 

We consider a single QD inside a semiconductor mi- 
crocavity that is pumped with a laser pulse and inter- 
acts with a cavity mode. It is assumed that the system 
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is initially prepared in its ground state. We consider a 
solid-state reservoir for the exciton population and we 
focus on time scales which phonon effects are important. 
We neglect other sources of damping in the system. We 
model the QD by a two-level system with ground state \g) 
(the semiconductor ground-state) and first excited state 
|e) (a single exciton), separated by an energy huj ex . The 
phonon environment is modelled by a bath of harmonic 
oscillators of frequencies Uk, with the wavevector k. The 
Hamiltonian of the total system in the rotating wave ap- 
proximation is written as 

H = huj ex a ee + huj c a) a + hwkbjjbh 

k 

+ hg(a eg a + a f a ge ) + hf(t)(a eg +a ge ) 

+ a ee J2^k{bk + b f k ), (1) 

k 

where <7jj = a (o^) and bk (b k ) are the annihilation 

(creation) operators for cavity mode, and fcth phonon 
mode, respectively. The parameter g is the coupling con- 
stant of the exciton and cavity mode, and f(t) is a real 
envelope function of the driving pulse. The last term 
in the Hamiltonian describes the exciton-phonon inter- 
action. In this term, is the corresponding coupling 
constant. The coupling of the confined exciton to the 
acoustic phonons by means of the deformation poten- 
tial tends to dominant the dephasing dynamics, over the 
piezoelectric interaction or coupling to optical phonons 
[19( 1 . In this case, the coupling constant is given by 
A fe = kD(k)y / 2nuj k V [HJ where n is the sample density 
and V is the unit cell volume. D(k) is the form factor of 
the confined electron and hole in the ground state of the 
QD. The Hamiltonian in the interaction picture can be 
written as 



Hint = H + Hj 



(2) 



where we decompose the coherent-field part and environ- 
ment part as follow 

H = hg(a eg ae lAt + al& ge e- lAt ) 

+ hf{t)(a eg e^ 1 + a ge e- iu ^), 
Hr = ZeeJ2 X k(he- iUht + Pe^). (3) 



In this equation A = u> ex — uj c is detuning between the 
exciton excitation energy in the QD and cavity field en- 
ergy. 

Now we consider the Liouville equation of density matrix 
in the interaction picture 



-^-j.[pt,H int \. 



(4) 



We define the reduced density matrix p for the exciton- 
photon system by tracing out the phonon degrees of free- 
dom in the total density matrix, p = Tr p h(pt)- Now we 



consider the master equation in the Born approximation 
[l|-[2| in the case of the phonon interaction while we con- 
sider the gain and pump parts exactly. Phonons are one 
of the slowest process and this kind of reservoir has a cor- 
relation time of the order of a few picosecond and this 
reservoir is naturally non-Markovian. To consider non- 
Markovian dynamics we have used time convolutionless 
projection operator method 20], up to second order of 
expansion. We assume an uncorrelated state for initial 
state of the exciton-photon system and phonon reservoir. 
At the initial time t = the phonon system is assumed to 
be in a thermal equilibrium at temperature T. Then the 
density operator of the exciton-photon system satisfies 
the following dynamical equation 



p(t) = -[ P (t),H ] 



{[a ee ,a ee p(t)]K(t-t') 



[a ee , p(t)a ee ]K* (t - t'))dt' . 



(•5) 



The first term describes the coherent evolution of the 
density matrix p under the action of the Hamiltonian H 
of the dot-cavity-pulse system. The kernel K which is 
the correlation function of the environment is written as 



1 f c 



du>j(uj) 



coth( 



2k B T 



)cos(u)t) — isin{ujt) 



(6) 

with Boltzmann constant ks- is the spectral density 
of the phonons which completely describes the interac- 
tion of exciton and phonons [22] . Here, we introduce the 
following spectral density 



i(^) = ^ A ^( w ~ Wfc )- 



(7) 



The density matrix dynamics is obtained under the Born- 
Markov approximation for exciton-phonon interaction 
and the strong exciton-photon interaction and pump ef- 
fects are described exactly. We can extract exciton dy- 
namics and photon statistics from this equation. 



III. EXCITON DYNAMICS UNDER A 
DRIVING PULSE 

In this section we consider the optical coherence of a 
driven QD under a pump pulse. Here we neglect the cav- 
ity mode and we consider optical coherence and exciton 
population dynamics under pulse excitation and effects 
of physical parameters such as pulse duration on these 
physical quantities. Then the density matrix of the exci- 
tonic system satisfies the following equation of motion 



Pex{t) = -[p ex {t),h(a eg a(t) + a ge a*(t))} (8) 



{[a ee ,a ee p{t)]K{t-t') 



[a ee ,p(t)a ee ]K*(t-t'))dt', 



(9) 
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where a(t) = fii)^^'* 1 . Exciton population and optical 
induced coherence in the QD system are defined through 
the different matrix elements of the density matrix. Ex- 
citon population and optical coherence are defined with 
the following set of equations, respectively 



P(t) = ia(t)(2N e (t) - 1) - P(t) / K(t - t')dt', 



N e (t) = 2iIm(a*(t)P(t)), 



(10) 



where P{t) = (e\p ex (t)\g) and N e (t) = (e\p ex (t)\e) . We 
assume at t = the QD be in its ground state and at 
this time it is excited with a Gaussian pulse excitation 

with envelop function f(t) = -^=-e~^ where a is the 
pulse width and A is a measure of pulse amplitude. For 
numerical integration of this set of equations, we shall 
take a GaAs QD with a spherical shape. In this case the 
spectral density is given by 



(cr e - a h f 
4ir 2 pc 5 



w 3 e 



(11) 



where a e and ah are the bulk deformation potential con- 
stants for electron and hole, c is the sound velocity in 
the sample and I is the electron and hole ground-state 
localization length (we assume a spherically symmetric 
harmonic confinement potential for the QD and electron 
and hole in the ground state). We use the following nu- 



merical values cr e — (j/j = 9eV, p = 5350 



kg 



5150^ 



and I = 4.5nm (these material parameters are approxi- 
mately acquired from [23[ ). Figure ^ shows plots of the 
time evolution of the exciton inversion for two values of 
pulse duration. In first picoseconds of dynamics the time 
evolution shows a strong decrease of exciton inversion due 
to the phonon effects and then we see a stable oscillation 
in inversion behavior during the pulse duration. It is 
clear from the figure that the phonon effects can prevent 
exciton generation. On the other hand, we see the com- 
plex behavior on the same timescales of initial dynamics 
for each pulse duration and after that small oscillations 
will continue at the end of pulse duration. Then we con- 
clude that in the first steps of dynamics the influence of 
phonons is a very important damping effect. Figure© 
shows plots of ImP(t) to consider the time evolution of 
optical coherence. As in the case of exciton population, 
optical coherence experiences a very rapid decrease dur- 
ing some first picoseconds. After this strong decrease we 
see a very small stable oscillations in optical coherence. 
Therefore, we concloud phonon effects are very impor- 
tant on timescales smaller than the spontaneous decay 
time and we can consider phonon reservoir as dominant 
damping source during the first steps of dynamics. 



IV. INTERACTION OF QD WITH CAVITY 
MODE 



case, the density matrix for the system satisfies Eq.©. 
By using Eq. ([5]) one can get a set of differential equations 
that describe the evolution of the populations and coher- 
ence of the cavity-QD system. In the basis of product 
states between the QD states and Fock states of the cav- 
ity mode (| en), \gn)) we calculate the matrix elements 
of the exciton-photon density matrix. By taking the ma- 
trix elements in Eq. ([5]) we get the following set of linear 
differential equations for the populations and coherence 
in the QD-photon system (we have used the notation 
Pin,jm — {i n \p\j m ) hi which i and j refer to QD states) 



Pen— l,en— 1 

(<) = igy/n(p en -i,gn(t)e 



-iAt 



if{t)(Pen-l,gn-l(t)e - Pgn-l.en-lit)^" 1 



Pgn,gn(t) = wVn{Pgn,en-l(t)e lAt - Pen-l,gn(t)e lAt ) 
+ if^){P 9 n,en{t)^ xt ~ Pen.gn^-^), 
Pen-l,gn(t) = igVn(pen-l,en-l(t)e l ' 



Pgn.en — 1 (^)^ ) 

Ju! ex t\ 

(12a) 
(12b) 



- Pgn,gn{t)e l f ) — p en -l,gn (t) / K(t — t')dt', 

JO 

Pen-l,gn-l{t) = {Pen-l,en-l [t)e™ ext 



(12c) 



Pgn-l.gn-li^e^) - p en -l, gn -l(t) / K{t - t')dt' 



(12d) 

In the absence of pulse excitation, the matrix elements 

Pen-l,en-l(t), Pgn,gn(t), Pen-l,gn(t) and p gn ,en-l(t) , for 

a given photon number, satisfy a closed set of differential 
equations. However, the excitation pulse couples the dif- 
ferent terms to each other and an infinite set of equations 
has to be solved. In the process of obtaining the above 
set of equations we neglect the terms like p gn ,gn-i(t) 
and p e n,en-i{t) because these terms do not have physical 
meaning related to the conditions under consideration. 
These terms show a coherence in photon field while the 
QD remains in its state. This could be related to photon 
damping which we have neglected such kind of terms. On 
the other hand, we maintain terms like p en ,gn(t) which 
describe coherence in QD system while photon number 
is constant. As is clear from (|12p these terms cam be 
generated during the dynamics by the pump pulse. 

As initial condition we take at t = the QD in 
its ground state and cavity field in the vacuum state 
PgO,go(0) = 1) and all other elements of the density ma- 
trix equal to zero. For the numerical integration, the set 
of equations can be truncated at a given value, which we 
take it equal to 90 (this value is choose with this condition 
that the results not change with increasing the number 
of equation). 

Photon statistics and material characteristics such as 
inversion population and optical coherence can be ob- 
tained from (fT2"j) . At first we consider Mandel parameter 
of the cavity field which is defined as [24| 



In this section we consider the interaction of the QD 
embedded in a microcavity with cavity mode. In this 



(13) 
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This parameter vanishes for the Poisson distribution, 
is positive for the super-Poisson distribution (photon 
bunching effect), and is negative for the sub-Poisson dis- 
tribution (photon anti-bunching effect). The mean num- 
ber of photons in the cavity is (other moments of n can 
be calculated in the same manner) 

(n) = ^ n [p en, en 

(«) + 

Pgn,gn 

(*)] ■ (14) 

n 

Mandel parameter for the case of resonant interaction 
(A = 0) and in the presence of detuning is plotted, re- 
spectively, in figures ([3]) and (® for two different values of 
pulse duration. As is seen, the cavity field mode exhibits 
non-classical (sub-Poissonian statistics) in the course of 
time evolution. Another important feature of this plot 
is the oscillatory behavior of Mandel parameter for time 
scales approximately two times of pulse duration. There- 
fore, the emitted photons to cavity mode by QD in the 
course of the excitation duration can be reabsorbed by 
QD and re-excite the QD then after the end of pulse du- 
ration we have oscillations in photon statistics. On the 
other hand, it is clear that with increasing the detuning 
feature the amplitude of oscillations in Mandel parame- 
ter decrease. 

Another important quantity in photon statistics is sec- 
ond order coherence function, g^(t,r) [J, [HI which is 
a two-time correlation function. Here we consider this 
quantity for the case of zero time delay, g^ 2 '(t,r = 0). 
This quantity can be used as an indication of the possi- 
ble coherence of the state of the photon system. For the 
single mode cavity field g( 2 >(t,T = 0) has the following 
definition 

gW(t,r = 0) = \ J 15 

J2n n ( n - jj [Pen.enjt) + p gn ,gn(t)] 
(E„ n \P 

In the case of resonant interaction and off-resonant inter- 
action, the plots of this quantity are shown in figures© 
and ([5]), respectively. The figures show non-classical na- 
ture of emitted photons (photon anti-bunching). This 
quantity shows similar oscillatory behavior to the Man- 
del parameter and its oscillatory behavior continue up to 
times twice the pulse duration. According to these plots 
the detuning effects on g^ 2 '(t, r = 0) are similar to its ef- 
fects on the Mandel parameter and cause the amplitude 
of oscillation be reduced. Therefore, in this conditions 
without any restriction on physical parameters (damping 
coefficients and coupling constant) it is possible that QD 
emits anti-bunched photons with sub-Poissonian statis- 
tics. The possibility of emitting anti-bunched photons 
with sub-Poissonian statistics by a single QD has been 
considered experimentally [25] . 

The time evolution of the QD coherence in the process 
of one photon interaction P(t) = (e0|p(t)|<?l) is shown in 
figures® and (O for different values of pulse duration 
and detuning. In these figures we plot imaginary part of 



P(t). These figures indicating occurrence of decoherence 
(damping of the imaginary part of polarization) in the 
system. The main source of this decoherence is phonon 
interaction. In the case of pulse with long duration we 
see an irregular oscillation in some time periods. It is 
clear that detuning prevents the coherence in this system. 
However the detuning is increased the imaginary part of 
coherence P(t) and increasing of detuning leads to the 
regular oscillatory behavior and causes damping will de- 
crease. In turn, because of the detuning, which weakens 
the dynamics, the pumping should be increased. Hence 
these two parameter can be considered as some experi- 
mental parameters for controling the decoherence in the 
QD systems on the timescales under consideration. On 
the other hand, by comparison Fig.® with Fig.® we 
can conclude that while the QD interacts with a cavity 
mode its optical coherence between energy levels has a 
longer life time. Then this can be considered as another 
experimental condition for controlling of optical coher- 
ence. 



V. CONCLUSION 

In this paper we have considered phonon effects (de- 
phasing effects) on optical properties of a pulse driven 
QD. We have shown that these effects strongly affect the 
Rabi oscillations and optical coherence. In the time scales 
which spontaneous emission and non-radiative recombi- 
nation do not play an important role in the dynamics 
(characteristic times of these effects are much longer than 
the characteristic time of phonon reservoir) the phonons 
strongly affect optical properties of QD. In the case of 
the interaction of system under consideration with cav- 
ity mode we have shown that emitted photons are anti- 
bunched and obey the sub-Poissonian statistics. Then 
in the microcavity with high quality factor which con- 
tains a single QD it is possible to generate non-classical 
light in the first some ten picoseconds. Here, we have 
considered a Gaussian pulse as a pump. We have shown 
that with the ending of pump, oscillations in the pho- 
ton statistics continue until times twice the pulse dura- 
tion. This relates to cavity photon which remains in the 
cavity and after ending of pump re-excites the QD. On 
the other hand, we have considered the detuning effect 
on the optical coherence of QD and we have seen that 
detuning can prevent decoherence effects. Hence, detun- 
ing can be considered as a controlling parameter of op- 
tical coherence. While QD interacts resonantly with the 
cavity mode, we have found that its optical coherence 
has a longer life time in comparison with its interaction 
with classical pulse. Then by putting the QD in the cav- 
ity it can maintain its coherence between energy levels. 
Therefore, the off-resonant interaction of a QD with cav- 
ity mode can be considered as an experimental tool for 
suppressing decoherence effects on the exciton. 
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FIG. 1: Plots of exciton inversion versus time for two different values of pulse duration: (a) a = lOps, (b) a = 40ps. Material 
parameters are pointed out in the text and T = 30K. 
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FIG. 3: Mandel parameter versus time for pulse duration a = lOps and T = 30K for two different values of detuning 
A = 0.0, A = 1.0. 
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FIG. 4: Mandel parameter versus time for pulse duration a = 40ps and T = 30K for two different values of detuning 
A = 0.0, A = 1.0. 
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FIG. 5: g (2) (t,T = 0) as a function of time for pulse duration a = Wps and T = 30A' for two different values of detuning 
A = 0.0, A = 1.0. 
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FIG. 6: g {2) (t,T = 0) as a function of time for pulse duration a = 40ps and T = 30A' for two different values of detuning 
A = 0.0, A = 1.0. 
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FIG. 7: ImP(t) as a function of time for three different values of detuning and for pulse duration a = lOps. In this plot 
T = 30K. 
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FIG. 8: ImP(i) as a function of time for three different values of detuning and for pulse duration a = 40ps. In this plot 
T = 30K. 



